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We present a shear-transformation-zone (STZ) theoretical analysis of molecular-dynamics sim- 
ulations of a rapidly sheared metallic glass. These simulations are especially revealing because, 
although they are limited to high strain rates, they span temperatures ranging from well below to 
well above the glass transition. With one important discrepancy, the simplified STZ theory used 
here reproduces the simulation data, including the way in which those data can be made to collapse 
approximately onto simple curves by a scaling transformation. The STZ analysis implies that the 
system's behavior at high strain rates is controlled primarily by effective-temperature thermody- 
namics, as opposed to system-specific details of the molecular interactions. The discrepancy between 
theory and simulations occurs at the lower strain rates for temperatures near the glass transition. 
We argue that this discrepancy can be resolved by the same multi-species generalization of STZ the- 
ory that has been proposed recently for understanding frequency-dependent viscoelastic responses, 
Stokes-Einstein violations, and stretched-exponential relaxation in equilibrated glassy materials. 
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I. INTRODUCTION 

A remarkable scaling property of viscosity as a func- 
tion of shear stress has been observed by by Olsson and 
Teitel [l| in simulations of a strongly overdamped, ather- 
mal, two-dimensional, amorphous system near its jam- 
ming transition. More recently, an approximately simi- 
lar scaling behavior has been seen by Guan, Chen and 
Egami 2] (GCE) in molecular dynamics simulations of 
a rapidly sheared, three dimensional metallic glass that 
undergoes a thermodynamic glass transition. These re- 
sults suggest that there may be a useful, minimal descrip- 
tion of the underlying physics for the seemingly complex 
dynamics of flow in glass- forming fluids, and that some 
scaling properties of a zero-temperature jamming transi- 
tion may persist in rapidly deforming systems far from 
jamming conditions. 

The distinguishing feature of the GCE simulations is 
that they cover an unusually wide range of system param- 
eters. They are carried out at kinetic temperatures rang- 
ing from 100 i^ to 1100 ii', i.e., from well below to well 
above an apparent glass transition. Importantly, they 
are carried out at strain rates from 10^ s^^ to 10^^ s~^, 
all well above values ordinarily accessible in the labora- 
tory. By driving their systems at large strain rates, GCE 
achieve steady-state, nonequilibrium shear flows at tem- 
peratures far below the glass transition. Under stress, 
molecular rearrangements become faster and the internal 
disorder increases, while the kinetic temperatures remain 
low. GCE are not probing truly glassy behavior in this 
way; however, they are discovering important features of 
glass dynamics by looking from this point of view. 

In this paper, we explore the implications of the GCE 
results by comparing them with the predictions of a 
shear-transformation-zone (STZ) theory ^3^. This com- 



parison is timely because Voigtmann [6| recently has pub- 
lished an analysis of the same data usingan extended 
version of mode-coupling theory (MCT) [S-Q- These two 
theories of amorphous plasticity would seem at first sight 
to have non-overlapping regions of validity. 

MCT starts with a liquidlike description of a many- 
body system. Its strength is that the coupling terms that 
emerge from its approximate closure of the many-body 
equations of motion can be evaluated in terms of observ- 
able structure factors, thus basing the theory directly 
upon microscopic dynamics. Its weakness, however, is 
that this closure approximation is accurate only if the 
density fluctuations are distributed Gaussianly. There- 
fore, MCT is effectively a mean-field theory in which 
motions are determined by averaging over large numbers 
of weakly correlated events. Such an approximation be- 
comes qualitatively incorrect at low temperatures, where 
the flow is governed by sporadic, thermally activated 
events that are rare but intense [Ifll- Voigtmann's results 
may imply that MCT retains some phenomenological va- 
lidity for flowing states in the low-temperature regime; 
but they cannot in principle have the same fundamental 
validity that MCT can claim at higher temperatures. 

In contrast, the STZ theory starts with a solidlikc pic- 
ture in which thermally activated, localized flow defects 
play the central role in controlling plastic deformation. 
It predicts that, at low temperatures, the yielding tran- 
sition as a function of stress is an exchange of dynamic 
stability between jammed and fiowing states; and it de- 
scribes the mechanical behavior of the system on both 
sides of that transition. When combined with an equa- 
tion of motion for the effective disorder temperature of 
the configurational (i.e. structural) degrees of freedom, 
the STZ theory successfully predicts shear-banding insta- 
bilities in agreement with simulations [ll|, II4I and, as will 



be seen here, describes a transition to liquidlike behavior 
at large stresses. It is in the latter regime, however, that 
the STZ theory encounters a weakness complementary 
to that of MCT. The highest GCE strain rates are com- 
parable to inverse molecular relaxation times, where the 
system must become a true fluid, and where the flow de- 
fects are no longer dilute and dynamically independent of 
each other. Nevertheless, like MCT at low temperatures, 
the STZ theory seems to retain at least a phenomenolog- 
ical validity up to quite high strain rates. Why does this 
happen? Is there a region of overlapping validity of the 
MCT and STZ theories? Is this theoretical agreement 
related to the approximate scaling behavior observed by 
GCE? 

Our tentative answer to these questions is that, to 
a first approximation, our results are sensitive to only 
a few features of the molecular interactions. The ob- 
served behavior is governed largely by general thermo- 
dynamic principles supplemented by dimensional analy- 
sis. The relation between stress and strain rate in the 
rapid-deformation limit is determined almost entirely by 
the density of STZ's, which is determined by the effec- 
tive disorder temperature. The details of specific mech- 
anisms, such as the thermal activation rate for internal 
STZ transitions, are almost irrelevant in most, but not 
all, of the regime of interest. 

Interestingly, the agreement between the STZ theory 
and the GCE simulations fails for the smaller strain rates 
at temperatures near the glass transition ~ for reasons 
that we think we understand. All of the analysis pre- 
sented here is based on a single-species STZ model, which 
we know to be incorrect for weakly perturbed, well equi- 
librated, glass forming systems at low temperatures. Un- 
der the latter conditions, a statistical argument predicts 
that STZ's with a wide range of internal transition rates 
become thermodynamically stable, and that the STZ's 
at the slow end of this range control the viscous response 
|13| . This distribution of transition rates accounts accu- 
rately for the very broad response peaks observed in oscil- 
latory viscoelastic measurements, and for related exper- 
imental observations such as Stokes-Einstein violations 
and stretched-exponential relaxation 



141. Indeed, the 



success of the STZ theory in predicting these intrinsically 
glassy, near-equilibrium phenomena gives us confidence 
in the present attempt to extend the analysis to strongly 
driven situations; but it also tells us where to expect dif- 
ficulties. 

We start, in Section |lll by presenting the steady-state 
STZ equations of motion in a more general and compact 
form than the one that has appeared in previous publica- 
tions j3ji5|]. Then, in Sec. lIIIl we summarize the primarily 
kinematic arguments that lead to our equation of motion 
for the effective temperature. In Sec. lIVi we describe the 
ways in which we use the simulation data to evaluate the 
theoretical parameters. Section |V] is devoted to our in- 
terpretation of the GCE scaling analysis. We conclude 
in Sec. I VII with a discussion of open questions. 



II. STZ DYNAMICS 

The basic premise of the STZ theory is that, in a closely 
packed, amorphous material, deformation is enabled by a 
fiuctuating population of rare, localized, two-state, flow 
defects, i.e. STZ's. Let n± denote the number densities 
of STZ's oriented parallel and antiparallel to a deviatoric 
shear stress s. These densities satisfy a master equation 
of the form 

To n± ~ R{±s) rizfi — i?(Ts) n± 



+ f(s,x) [noo(x) -n±\ 



(2.1) 



On the left-hand side, the factor tq is a microscopic time 
scale, of the order of picoseconds for metallic glasses. On 
the right-hand side, the first terms containing the factors 
R{±s) are the rates of forward and backward STZ tran- 
sitions (in units oi Tq ). In general, these rate factors 
are functions of both the ordinary and effective temper- 
atures. They determine the plastic strain rate 7^', which 
we write in the dimensionless form: 
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(2.2) 



where vq is the average volume per molecule, and eg is a 
dimensionless constant of the order of unity. 

The second terms on the right-hand side of Eq. (|2.ip 
are the rates of STZ creation and annihilation. These are 
fluctuation- activated processes, expressed in the form of 
a detailed-balance relation in which the steady-state STZ 
density is determined by the effective temperature %: 
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(2.3) 



Here, x is written in units of the STZ formation energy 
ez, i.e. X — ks Teff/ez- To a first approximation, ez is 
the energy required to create the amorphous analog of a 
vacancy-interstitial pair. Note that the effective temper- 
ature introduced here is the one that is defined thermo- 
dynamically in (15| . We presume that it is the same as 
the effective temperature determined by a fiuctuation- 
dissipation relation, as described in [16]] and references 
cited there; but we know of no rigorous proof of that 
assertion. 

T{s,x)/to is an attempt frequency consisting of addi- 
tive thermal and mechanical parts: 



r(s,x) = p + r(s,x). 



(2.4) 



The term p = p{0), where 9 = ksT/ez, is best under- 
stood as a dimensionless, thermal noise strength. It is a 
super- Arrhenius factor whose strong temperature depen- 
dence governs the equilibrium glass transition, and for 
which one of us ( JSL) has proposed the outlines of a first- 
principles derivation. (iTl [l8j We expect that p{9) is of the 
order of unity at high temperatures, and that it decreases 
rapidly toward zero as 9 falls through the glass tempera- 
ture. In the absence of external driving, p(9) controls the 



rate at which the system undergoes structural aging. In 
much of the hterature, the quantity p{d) exp [—1/0) /tq 
is defined to be the a relaxation rate r~^. However, as 
pointed out in [14!|, this definition is not the same as the 
conventional assumption that Tq, is directly proportional 
to the viscosity. 

In analogy to p{0)/to, the quantity r(s, x)/to is the 
contribution to the attempt frequency in Eq. (j2.4p due to 
mechanically generated noise. As in earlier papers |5l.ll9l. 
l2q |. we assume that T{s, x)/to is proportional to the rate 
of entropy production. In steady fiow, all of the work 
done on the system is dissipated as heat; therefore, the 
rate of energy dissipation per unit volume is 2 7^' s. To 
convert this rate to a noise strength with dimensions of 
inverse time, we multiply by the volume per noise source, 
i.e. the volume per STZ, vq exp(l/x), and divide by an 
energy conveniently written in the form 2 ep sq vq . The 
fact that the stress sq, introduced here for dimensional 
reasons, turns out to be the low-temperature yield stress 
has been one of the more interesting surprises in this 
theory. The resulting formula for the noise strength is 



Then, knowing r(s, x), we can use Ea. (l2.5p to write: 
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(2.5) 



The GCE simulations involve only steady-state defor- 
mations; thus all of the information that we need is con- 
tained in the stationary solutions of Ea. (|2.1[) . obtained 
by setting n± = and solving for the n±. The resulting 
expression for the strain rate in Ea. (l2.2[) is 



q = eoT{s,x)e 
where 
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r(s,x), (2.6) 
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and 
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(2.7) 



(2.8) 



According to Eq. (|2.8l) . T{s,x) is the bias between for- 
ward and backward transitions. The second law of ther- 
modynamics (see [l5|) requires that 



T{s, x) = tanh 



vgs 
ezX 



(2.9) 



With Eq. (|2.5p . Eq. (|2.6p becomes a quadratic equation 
that can be solved for q or, more conveniently, for r(s, x)- 
The result is 

r(s,x) = I Q{s,x) + \ VQG^,x)' + 4p(0)Q(s,x), 

(2.10) 
where 

Q(5,x) = 2C(5)r(s,x) - -2C{s) - p{e). (2.11) 
so 



^e-Vxr(,,;,). 



(2.12) 



Note that r(s, x) is a non-negative, symmetric function 
of s that vanishes like s^ at s = 0. 

In the athermal limit where p = 0, both T{s, x) f^nd q 
vanish for s < Sy (the jammed state); whereas, for s > Sy 
(the fiowing state). 



g = 26oe-i/^C(s)[r(s,x)-- 
L s 

The yield stress Sy is the solution of 

SyT{Sy,x) = So- 



(2.13) 



(2.14) 



Therefore, Sy ^ sq when Sy is large and T(sy,x) — 1- 
For small but nonzero p, according to Ea. (j2.10p . both 
r(s,x) and q make smooth transitions near s — Sy be- 
tween viscous and flowing states. 

We also need to examine the viscous limit at small q, 
where r(s,x) <?C p and x ~ ^- Keeping only the linear 
term in the relation between q and s, we find 

—'»'-'" (^cwfk) 51- <-^' 

The Newtonian linear viscosity, for vanishingly small q, 
and expressed here in units of stress, is r/jv = s/q. Note 
that this formula produces the conventional result in 
which rjN ex Tq multiplied by a slowly varying function of 
9. It is this result that, according to [ij], substantially 
underestimates t]]^ at low temperatures and, therefore, is 
responsible for the important discrepancy between this 
theory and the GCE simulations. 

C{s) is necessarily a symmetric function of s. Wc write 
it in the form: 



C{s) — exp 



_^g-.V24 



1/4 



(2.16) 



On the right-hand side, the first factor in square brackets 
is a thermally activated (Eyring-like) rate in which a bar- 
rier of height 9e is reduced quadratically by the applied 
stress via a Gaussian factor, thus preserving the required 
symmetry. In what follows, we find that we can approxi- 
mate this term by unity. The second factor is similar to, 
but not strictly the same as a Bagnold rate, proportional 
to the square root of the stress at large s for dimensional 
reasons. It does seem to be significant at low tempera- 
tures and high strain rates. The only quantity that can 
play the role of sb in this formula is the pressure, which 
may be roughly proportional to the temperature. Thus, 
we write sb = (T/Tb) GPa, where the "Bagnold tem- 
perature" Tb is a constant. 



III. EFFECTIVE-TEMPERATURE 
KINEMATICS 

Turn now to the effective temperature x- Iii general, 
X is the temperature of the configurational degrees of 
freedom of the system, thus it characterizes the system's 
state of structural disorder. Its role here is primarily to 
determine the density of STZ's, as indicated by its ap- 
pearance in the Boltzmann factor in Eq. (|2.3p . We pro- 
pose that X is determined mostly by kinematics rather 
than by specific details of the molecular interactions, and 
that this property of x is largely responsible for the uni- 
versality observed at high strain rates. 

At low temperatures, where changes in the state of 
glassy configurational disorder can be induced only by 
externally driven deformation and not by thermal fiuc- 
tuations, i.e. where p = 0, there must be a direct rela- 
tion between the dimensionless strain rate q and the low- 
temperature effective temperature, say x = xil)- When 
the strain rate is much smaller than any relevant relax- 
ation rate in the system, so that (7 <C 1 , then dimensional 
analysis based on the fact that there is no intrinsic rate 
comparable to (//tq requires that xi^) be equal to a con- 
stant, say xo- 

In the opposite limit, when q is of the order of unity, 
we assume that the relation between q and x has an Ar- 
rhenius form, q ^^ qo exp(— ^/x). Both x and ^ are 
energies measured in units of ez- If ez is the only en- 
ergy scale in the system relevant to configurational re- 
arrangements, then we expect that A is of the order of 
unity. Similarly, if Tq is the only intrinsic rate, then 
go also should be of the order of unity; but here there is 
an additional uncertainty about how accurately we have 
estimated tq. We return to these estimates in Sec. IIVI 

The relation between strain rate and configurational 
disorder was discovered in numerical simulations by Hax- 
ton and Liu 2^. It was discussed further in [2lj, where 
it was written in the form of a glasslike relation between 
a "viscosity" q~^ and the temperature x: 
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and, in analogy to Vogel-Fulcher-Tamann (VFT), 



ae//(x) 



Xi 



X-Xo 



exp 



X-Xo 



XA -Xo 



(3.2) 



Thus, X ~^ Xo in the limit of small strain rates, and 
X ^ cx) as g — >■ go- The exponential cutoff in Eq. p.2p 
is needed in order that the VFT divergence at small x 
transform smoothly to the Arrhenius law at large x- Pre- 
vious calculations have used 6 = 3. Experience with these 
formulas, as in [2l|, leads us to conclude that they are 
more cleanly reliable than their VFT analog for viscosity 
as a function of ordinary temperature. (In this connec- 
tion, see the calculation of the viscosity in [14|.) 

The equation of motion for x itself is a statement of 
the first law of thermodynamics; it describes entropy flow 



through the slow, configurational degrees of freedom into 
the fast thermal motions of the molecules. Near steady 
state, it has the form 



-i/x 



r(5,x) x('z)-x +'«p(^)(^-x) • (3.3) 



X oc e 



The first term in square brackets on the right-hand side is 
the rate at which x is driven toward x('?) by the mechani- 
cal noise strength r(s, x)- The second term, proportional 
to p(0), is the rate at which thermal fluctuations drive x 
toward the ambient temperature B. k is a dimensionless 
parameter of the order of unity. Since x is a measure of 
conflgurational disorder, we see here explicitly how p{ff) 
controls the rate of structural aging in undriven systems. 
The competition between these two terms in Eq. p.3p de- 
termines the value of x; it is close to x(g) for large r(s, x), 
and close to Q when the system is driven slowly so that 
r(s,x) is small. 

There is one complication that must be dealt with at 
this point. Equation (j3.3|) . as written, implies that the 
steady-state x niust lie in the interval between x('Z) and 
Q. If we assume that x('Z) ~ Xo is a constant for small 
enough g, then a system initially prepared with > Xo 
would be "cooled" to x < ^ when driven at a small strain 
rate. This behavior seems implausible; so far as we know, 
it is not seen in other simulations, e.g. [26|, in which x 
is measured directly. In [2l| , this problem was corrected 
by setting Xo = ^ when Q exceeds xo, and by rescaling 
Xi and XA accordingly. Speciflcally, 



Xo 



Xo, 



for 9 < Xo, 
for 6* > Xo 



Xi 



for 9 < Xo, 
(xo/xo)xi, for 6* > Xo 



Xi, 



and 



^ XA, for 9 < Xo, 

^^ l(Xo/xo)XA, for 61 >xo 



In summary, we use Eq. (|3.3|) in the form 

Tis,x)xiq) + ^p{d)0 



X 



r(s,x) + Kp(0) 



(3.4) 
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along with Ea. (l3.ip . to determine one relation between 
X, g and s. We then use Eq. (|2.12p to compute both g 
and X as functions of s. 



IV. EVALUATION OF PARAMETERS AND 
COMPARISON WITH THE GCE DATA 

In Fig. [1] we show a selection from the raw 
data of GCE. The data sets shown here are char- 
acteristic of the eleven such sets that we have used 
in these analyses. They are stress versus strain- 
rate curves for seven different temperatures, T = 
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FIG. 1: (Color online) Log-log plot of stress s versus di- 
mensionless plastic strain rate q — ro7^'. The top three 
solid black theoretical curves and the associated open-square 
data points are for temperatures T = 100 K, 300 K, and 
500 K. The middle dashed blue curves are for T = 700 K 
(blue circles), and 840 -ft' (blue triangles). The bottom two 
solid red theoretical curves and open-circle data points are 
for T — 94:0 K and 1100 -ft". Agreement between theory and 
simulation is good for the low temperatures (black curves) 
at the top and for the high temperatures (red curves) at the 
bottom; but it fails at small strain rates for the intermediate 
temperatures near the glass transition. 



100 K, 300 K, 500 K, 700 K, 840 K, 940 K and 
1100 if, reading from top to bottom. The data points 
are shown here along with our best-fit theoretical curves. 
The GCE time constant is tq = 0.1 picoseconds; there- 
fore, the dimensional strain rate is 7^' = 10^^ q s~^. We 
have arbitrarily chosen the upper limit of validity for the 
STZ theory to be at 9 = 10"^, i.e. at 7P' = 10" s^\ Be- 
yond this point, the simulations show strain-rate soften- 
ing - possibly an indication that the system is liquifying 
and is no longer consistent with the solidlike STZ model. 
Our strategy for choosing parameters has been to look 
for the simplest possible version of the theory that is 
consistent with the data. Accordingly, our starting as- 
sumption is that the Eyring-like rate factor is negligi- 
ble; that is, we set 9e = in Eq. (|2.16p . We have tried 
various nonzero values of 9e , and have found no overall 
improvement in the results. At low temperatures, only 
large stresses - above the yield stress - come into play, 
reducing the activation barrier via the Gaussian factor 
in Eg. (12. 161) . Conversely, at high temperatures, the driv- 
ing stresses may be small, but the activation barriers are 
small compared to /csT. The exceptions are at the inter- 
mediate temperatures, for small stresses and small strain 
rates, where the theory breaks down for the more inter- 
esting reasons mentioned in Sec. ID If we choose 9e large 



enough that, just by itself, it substantially increases the 
viscosity and improves the fit, for example, at T = 840 K 
(the lower dashed curve in Fig.[T]), then we qualitatively 
ruin the fit at T = HOOK. Ultimately, the Eyring bar- 
rier must be important at smaller strain rates and low 
temperatures; but we see no need for it here. 

The dimensionless constants eo in Ea. (l2.2p and n in 
Eq. p.3p should be of the order of unity. In the absence 
of better information, we set eo = k = 1 . Similarly, we 
set ez/vQ = 1 GPa and, in the next paragraph, use the 
high-temperature viscosity to check that this estimate is 
reasonable. 

At high temperatures and small stresses, with C{s) ~ 
^oVq/cz ~ 1, Equation (|2.15l) tells us that the linear 
Newtonian viscosity is 



m = [^ + pm 



2p{e) 
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(4.1) 



The last approximation is valid for temperatures high 
enough that p ~ 1. At T = 1100 if (the bottom curve 
in Fig.p, the GCE resuh is riN = 150. Then, Eq.(|33]) 
implies that 9 = 0.15, and therefore that Tz = ezj^B ~ 
lOOOK. This means that v^/ez ~ 10 uq, with uq mea- 
sured in cubic nanometers. Thus, a length scale of about 
half a nanometer gives us ez /vq « 1 GPa as in the pre- 
ceding paragraph. 

Now consider low-temperature situations for which the 
GCE data appear to indicate well defined yield stresses 
so(^)j implying that p{9) = to within the accuracy 
of the measurements. The temperatures for which this 
is true are T = \00K, 300 if, 500 isT and 600 if. In 
this regime, T{s,x) ~ 1; thus, in the limit of small q, 
Eq. (P?T^ becomes 
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(4.2) 



which we have used to make a first estimate of xo- We 
then have found that we can fit all four of these data sets, 
for the whole range of strain rates 10^^ < q < 10^^, 
by using the observed values of so(0), and by choos- 
ing Xo = 0.1, xi ^ XA = 0.2, A = 1.3, (Jo = 5, and 
6 = 3 in Egs. p.ip and p.2p . We also have chosen the 
Bagnold temperature defined following Eq. (|2.16l) to be 
Tb = lOOif. With this value of Ts, the Bagnold term is 
negligible for all but the lowest temperatures and largest 
stresses; but it does seem to be relevant in that regime. 
We have no reason to believe that this set of temperature- 
independent parameters is unique or optimal. 

The only remaining parameters to be chosen are the 
quantities p{9) and sq{9) for temperatures near and 
above the glass transition. Rather than trying to pre- 
dict these quantities theoretically, we have "measured" 
them by fitting the GCE data. The results are shown 
below in FigsI3]and[5l and are discussed in Sec. IVl 

The theoretical fits to the data in Fig. [T] are accurate 
for both the lower and higher temperatures, indicated by 
solid black and red curves respectively. However, there 
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FIG. 2: (Color online) Log-linear plot of the theoretical ef- 
fective temperature x as a function the dimensionless strain 
rate q for temperatures T = 100 K through 600 K (dark black 
line), 840 ii' (dashed blue line), and 940 A' and 1100 /C (thin 
red lines, from bottom to top). 



are significant discrepancies at small strain rates for the 
intermediate temperatures T — 700 K and 840 K, shown 
by the dashed blue lines. (There is a similar discrepancy 
for T = 800 K, not shown in the figure.) At T = 840 K, 
for example, the theoretical stress starts to fall below the 
observed values when q becomes less than about 10"'*; 
and the system switches over to linear viscosity at smaller 
strain rates. The multi-species STZ theory proposed in 
\13i [l4| predicts that the same thing will happen when- 
ever p is nonzero, but that the transition will occur at a 
smaller q and, therefore, at a larger viscosity. 

In Fig. [21 we show four examples of our theoretical x's 
as functions of strain rate. All four converge to the same 
large-g behavior shown in Eg. p. II) . The heavy black 
curve is the low-temperature function, x — x{q), valid for 
all cases in the range 100 K < T < 600 K where p ~ 0. 
It decreases logarithmically at small q toward xo = 0.1. 
The two red curves, for T = MO K and 1100 JiT, where 
p is large, fall quickly to their corresponding values of 
X — (^ = T/Tz at small q. The blue dashed curve is for 
T = 840 X, where p = 0.1. This curve falls substantially 
below the p — curve before levelling off at T/Tz — 0.12, 
again an indication that the small-q behavior is incorrect. 
The theory predicts too small a viscosity in this case; 
therefore, it predicts too small a value of r(s, x) relative 
to p in Eq. (|3.7l) . and thus predicts too fast an approach 
to the thermal limit. Note that the largest values of x 
shown in this figure are still small enough that the STZ 
density ~ exp (— 1/x) remains small. 
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FIG. 3: (Color online) Thermal noise strength p (open red 
circles) as a function of temperature T. 



V. DATA ANALYSIS 

In principle, an ideal glass transition occurs in an un- 
driven, equilibrated, amorphous system at a Kauzmann 
temperature Tq, below which the viscosity is infinite and, 
above which, the yield stress vanishes. The GCE scaling 
analysis implies that something like this transition is oc- 
curring in this rapidly deforming system. 

To begin to understand this situation, we show in Fig. [3] 
our "measured" values of the thermal noise strength p for 
the eleven different temperatures included in the GCE 
data. This function behaves more or less as expected, 
falling from a value of the order of unity at large T to 
become unmeasurably small below the glass transition. 
The most notable feature of this curve is the smooth 
transition to zero in the region 600 K < T < 840 K. 
For these temperatures, we have optimized the fits to be 
accurate at the higher strain rates, as shown in Fig. [1] by 
the dashed curves for T = 700 K and 840 K. The single- 
species STZ theory should be valid at large q. Therefore, 
these measurements of p must be taken seriously despite 
the fact that they appear to be inconsistent with VFT 
formulas for p of the kind developed in [13, [3 ■ 

The glass transition emerges more clearly if, instead of 
plotting p itself, we use the p-dependent formula for tjn 
in Eq. (|4.1[) to plot the inverse Newtonian viscosity 77^^ 
as a function of temperature. This is the function that 
GCE suggest can be fit by 
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otherwise 



(5.1) 



where Tq — 860 if and a — 1.23. Our result and the 
comparison with the GCE fit are shown in Fig. |4l Here, 
the fact that we are underestimating the viscosities for 
temperatures less than Tq does not make much difference 
because r;^ is already very small in that region. 

According to the STZ theory, there is no yield stress 
in the limit of zero strain rate when p is nonzero; the sys- 
tem must have a finite linear viscosity at small enough 
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FIG. 4: (Color online) Inverse Newtonian viscosity rj'^^ (open 
red circles) as a function of temperature T. Also shown (blue 
solid line) is the GCE fit to this curve given in Ea. H5.1|) . 

stresses and strain rates. Nevertheless, the theory pre- 
dicts that, at very small p, the system will appear to have 
a yield stress so long as the strain rate is not too small. 
Therefore, we can plot this apparent yield stress by plot- 
ting the observed stress at the smallest observed strain 
rates, as a function of temperature. The GCE data for 
q = 10~^ is shown in Fig. [SJ along with the GCE fit to 
this apparent yielding curve: 
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where f3 = 0.6. The yield stress predicted by the STZ 
theory would be given by the blue triangles in this figure 
for temperatures up to about 600 A' (for which p = 0), 
and then would drop abruptly to zero above that point. 
It would be interesting to see whether something closer 
to this behavior would appear if the GCE data were ex- 
tended to smaller strain rates. We also show our mea- 
sured values of the stress Sq in Fig- El In accord with 
the preceding remarks, Sq coincides with the yield stress 
at temperatures up to 600 K. At higher temperatures, 
as found in [2l[, sq levels off at about 0.55 GPa. This 
simple behavior is determined by the large-g dependence 
of the high temperature curves in Fig. [I] where the onset 
of nonlinearity is sensitive to sq- 

The power-law approximations shown in Eqs. (|5.ip 
and (I5.2p are similar to those used by Olssen and Tei- 
tel [ij in their analysis of a zero-temperature, jamming- 
unjamming transition controlled by the volume fraction 
of a model granular material. The GCE analog for the 
temperature controlled situation is 
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where r](s) ~ s/q (not the small-g Newtonian viscosity). 
A graph of 77^ as a function of st produced the data 
collapse shown in GCE Fig. 2. Our STZ version of that 
figure is shown here in Fig. [BJ where we have included 
the data for all eleven of the GCE temperatures. 
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FIG. 5: (Color online) Apparent yield stress from GCE data 
(open red circles) and the GCE fit to this data (dashed black 
curve) given in Ea. (|5.2p . The open blue triangles show the 
measured values of the stress sq. 



Note that the high-temperature (upper) branch in 
Fig. [S] has become nearly a single horizontal line at the 
constant inverse viscosity rj^ given in Eq. (j5.3l) . These 
lines do not lie exactly on top of each other in our graph 
because, as seen in Figs. |3] and [5l the power-law fits in 
Eqs. (|5.1|) and (|5.2p are quantitatively inaccurate. The 
low-temperature (lower) branch is not quite so clean, be- 
cause 77^ is small but nonzero for temperatures near Tq 
as shown in Fig. 2) The resulting failure of the scaling 
analysis appears in the three dashed curves in Fig. [B] for 
the temperatures T = 700 K, 800 K and 840 K, where 
the low-temperature behavior crosses over to an apparent 
linear viscosity with decreasing strain rate. This effect 
would be even more pronounced if, as predicted, the ac- 
tual yield stress drops to zero above 600 i^. Apparently, 
the data collapse seen in GCE does not indicate an exact 
scaling relation or the presence of some diverging length 
scale. Nevertheless, it remains an interesting qualitative 
feature of the data. 



VI. CONCLUDING REMARKS 

The analysis presented here makes it appear that 
the STZ theory of amorphous plasticity is realistic up 
to strain rates of about 10~^/ro, i.e. up to about 
one percent of the underlying molecular relaxation rate. 
Both the GCE simulations and our theoretical results 
exhibit an apparently simple scaling behavior across 
three decades of strain rates below that limit. Our 
STZ analysis implies, at least in a first approximation, 
that effective-temperature thermodynamics, rather than 
material-specific molecular interactions, control the sys- 
tem dynamics at high strain rates. 

As discussed in Sec.|Tl these STZ-based results are tech- 
nically beyond the range of the mode-coupling theory 
LD-Q, because they extend down to temperatures well 
below the glass transition where plastic flow is governed 
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FIG. 6: (Color online) Scaled inverse viscosity 77^ as a func- 
tion of the scaled stress st, computed from the STZ the- 
ory using all eleven GCE temperatures. The color scheme is 
analogous to that used in Fig. [l] The four solid red curves 
that collapse approximately to a single scaled inverse viscos- 
ity are for high temperatures in the range 900 AT — 1100 AT. 
The four solid black curves that collapse approximately to a 
single scaled yield stress are for temperatures in the range 
100 K — 600 K. The three dashed blue curves, which deviate 
from the scaling pattern, are for temperatures 700 K, 800 K 
and 840 K, reading from bottom to top. 



nonperturbative activation mechanisnis. However, see 
for a counter argument. Similarly, we do not see how 
our results can be reproduced in the framework of soft 
glassy rheology (SGR) ^ ^, because the STZ and 
SGR theories are so different kinematically, and because 
our present analysis is based so strongly on effective- 
temperature dynamics. Again, however, we point to [25|, 
which describes progress in reformulating SGR in a ther- 
modynamic framework analogous to that of [l5|. It will 
be interesting to see whether the results presented here 
do, in fact, distinguish the STZ theory from either of 
these two other competing points of view. 

Our analysis in this paper is based on several, fun- 
damental hypotheses that are potentially falsifiable, and 
which require further attention. First, in order to check 
the thermodynamic hypothesis, we need to test our the- 
oretical estimate of the effective temperature x by mea- 
suring it directly. Computational measurements of y in 
[221 served as the basis for the theory developed in [2l|, 
which led to the equations of motion for x presented here 
in Sec. IIIIl Analogous measurements of an effective tem- 
perature for a thermalized har d-sp here system have been 
reported recently by Haxton [2al; and an STZ analysis 
of Haxton's results has revealed interesting insights [23] . 



So far, we have not had a similar consistency check for 
our analysis of the GCE data. 

Second, our qualitative explanation of the discrepan- 
cies between the GCE data and single-species STZ theory 
is based on the multi-species reinterpretation of glassy 
viscosity proposed in [lj|. We need now to develop 
this reinterpretation into a quantitative theory, and test 
whether that theory corrects the discrepancies. This is 
our most unconventional effort. We have used large- 
strain-rate data, instead of the conventional small-strain- 
rate viscosity, to evaluate the thermal noise strength p{9). 
This procedure produces unexpectedly large values of 
p{6) at temperatures near Tq, as seen in Fig. |31 The 
argument in [ij| tells us that, because of the statisti- 
cally inevitable presence of "slow" STZ's, the viscosity 
increases more rapidly with decreasing temperature than 
would be predicted simply by the 9 dependence of p. 
This argument, when applied to undriven systems near 
their glass temperatures, accounts for Stokes-Einstein vi- 
olations and even for the stretched-exponential nature 
of various relaxation functions. The question is: Does 
it also account quantitatively for what look like anoma- 
lously large viscosities in the GCE data? 

A third class of questions arises naturally within the 
context of the first two. Our analysis indicates that the 
model of a metallic glass used in the GCE simulations un- 
dergoes some kind of glass transition at Tq = 860 K. As 
T decreases below Tq, the configurational noise strength 
p decreases smoothly and becomes unmeasurably small 
below about 600 K. As this happens, the system changes 
from a viscous, glass- forming liquid to a solidlike glass 
with an apparently well defined yield stress. At nonzero, 
steady, shear rates, this transition becomes smoother; the 
liquidlike behavior persists down to arbitrarily low tem- 
peratures with continuous increase in the stress level. 

Finally, an obvious question from a first-principles 
point of view is how the STZ transitions appear in the 
form of changes in the underlying atomic structure. Ref- 
erence [2^ , a sequel to the GCE simulation [3 , is a recent 
example of work along these lines. It is shown there that, 
during steady-state flow, the macroscopic Maxwell relax- 
ation time rj/Grx, (where Goo is the instantaneous shear 
modulus) is essentially the same as the time required for 
a localized topological change to occur in which an atom 
gains and/or loses one of its nearest neighbors. We ex- 
pect that further investigations of this kind will lead to a 
clearer microscopic picture of the STZ transition mecha- 
nism. 
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